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We present new results concerning the existence of static, electrically charged, 
, perfect fluid spheres that have a regular interior and are arbitrarily close to 

CN \ a maximally charged black-hole state. These configurations are described 

! | by exact solutions of Einstein's field equations. A family of these solutions 

■ had already be found (de Felice et al., 1995) but here we generalize that 

result to cases with different charge distribution within the spheres and show, 
<^ , in an appropriate parameter space, that the set of such physically reasonable 

solutions has a non zero measure. We also perform a perturbation analysis and 
| identify the solutions which are stable against adiabatic radial perturbations. 

We then suggest that the stable configurations can be considered as classic 
models of charged particles. Finally our results are used to show that a 
i . conjecture of Kristiansson et al. (1998) is incorrect. 
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H I. INTRODUCTION 

It is well known that static, spherically symmetric, uncharged perfect fluids cannot be 
held in equilibrium below a certain radius without developing singularities inside. This 
radius is 9M/4 for an incompressible fluid sphere of total mass energy M; it can be even 
larger than that for more realistic equations of state (Buchdahl, 1959). The possibility of 
holding a non-singular object in stable equilibrium but compact enough to be close (in fact 
arbitrarily close) to a black-hole state, is of great interest not only in order to judge the 
state of matter in this quasi-critical condition, that is being about to turn into a black-hole, 
but also to yield a classic model of charged massive particles which might have astrophysical 
and cosmological implications. Although one can reach this goal with non-perfect fluids, a 
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perfect fluid solution of the type mentioned was recently found (de Felice et al., 1995) but 
with the presence of an electric charge. Electric charges inhibit the growth of space-time 
curvature and therefore they are an efficient means of avoiding singularities inside matter 
and enhance stability even for quasi-critical configurations. 

In this paper we extend the above mentioned analysis (de Felice et al., 1995) to a more 
general case and, in Section 2, show the existence of regular solutions very close to a black- 
hole state as points in a parameter space (figures |T]-[3[). In Section 3 we study the stability 
of these solutions against adiabatic radial perturbations and show the domain of stable so- 
lutions in the same parameter space (figures |T|, |2|). Having solutions which describe very 
compact sources, we are able to verify the correctness of a conjecture (Kristiansson et al., 
1998) according to which the limit of regular embedding in the Euclidean space of confor- 
mally reduced space-like sections of the Reissner-Nordstrom metric, coincides with the limit 
of regularity of the internal solutions. In section 4 we show that this is not true. 

Throughout the paper we use geometrized units (c = 1 = G) and metric signature +2. 



II. EXISTENCE OF REGULAR QUASI-CRITICAL SPHERES 

In Schwarszchild coordinates the line element of a spherically symmetric space-time reads: 
ds 2 = -eMt 2 + e A dr 2 + r W + r 2 sin 2 Odd) 2 



where r\ = r](r, t) and A = A(r, t). The energy-momentum tensor of a perfect, electrically 
charged fluid takes the form (de Felice et al.,1990): 

Ti» = (p + p)u»u v + p<T + ^{F^F\ - -a^F aP F^) 

47T 4 

where F a ^ satisfies Maxwell's equations: 

F [a p M = 0, F% = **f. 

Here j M is the current four- vector, p and p are the energy density and the isotropic pressure 
of matter measured in its rest frame respectively; p and p are expressed in (geometrized) 
units of length^ 2 . 

In the spherically symmetric case the only non vanishing component of Maxwell's tensor 
is (Bekenstein, 1971): 

f t r = e _HsQ(r,t) 



r 2 



where: 

Q(r,t) = r e ^47rr 2 j t dr (1) 

<J 

is the total electric charge within a sphere of radius r at time t; it then follows that: 



F *f3 F 2( 1 2 
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where charge is expressed in units of length. If we require that the configuration be static, 
then = e - ^ 2 ^, hence Einstein's equations which one needs to solve for the interior 
metric, take the form: 

x fX 1\ 1 Q 2 
6 {-7-72)+^ = ^ + ^ (2) 

- Q ~ x {j + v) + h = %-^ p (3) 

where (') means derivative with respect to r. 

Let the sphere be incompressible with a constant matter energy density p m and the 
charge distribution within the sphere be given by: 

Q(r) = Qo (-)" (4) 

where n is a constant parameter, R is the coordinate radius of the boundary and Qo is the 
total charge. Integrating equation (|2]) with respect to coordinate radius r, one obtains 

6 =1 -T Pmr - (2n-l)R* - (5) 
Combining this result with (|l|) and we see that the charge density is given by: 



v t Qonr 11 - 3 ^ 0llr V 1 ~ T^ 1 - (2n-l)R2 



47rR n et 47rR n 

then we have n > 3 which prevents the divergence of charge density at the center. 
Introducing scale parameters such as: 

1 / 3 y /2 

po = 3 p - ro= (i^:J (6) 

and setting: 

Y=S-, { = - (7) 
Po r 

the Oppenheimer-Wolkoff equation and equation @ for such a sphere read respectively (de 
Felice et al., 1995): 

e -A = 1 _ 2 f - [a/(2n - l)]^^ (9) 

where 

An obvious condition is a > 0. Now we see that the set of equations (|3|), (§) and (|9]) is 
complete. 
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FIGURES 




FIG. 1. Configurations in the parameter space for n=5 



The boundary of the physical configuration is identified by the conditions: 

Y&) = 0, e^(6)=e- A (6) = l-2^-^ T ^ (n ~ 1) . (10) 

where |& = R/r Q . Then, for a given set of (a, n), one can solve equations (|3|), (H) and 
(§) and obtain the structure of the corresponding configuration (Note: When one integrates 
equation (|j) numerically from its boundary inward, pressure Y may diverge before £ reaches 
zero for some values of (a, £&, n). Since there is a singularity in such a configuration, 
it is physically meaningless. But we still use the parameter (a, n) to represent such 
solutions). But not all of the solutions are acceptable, two basic requirements of a regular 
static configuration are: 

1: That there be no singularity inside the sphere, 

2: That the radius of the boundary be larger than the external horizon size of the black 
hole with the same mass and charge. 




FIG. 2. Configurations in the parameter space for n=3 
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FIG. 3. Enlargement of figure Q around the meeting point 

For a given n, each point in the (£^, a)-plane represents a solution. In figures [L]-|3| we 
plot the functions which allow one to single out the points which correspond to physical 
configurations. 

Let M be the total mass energy of the solution, then the relations between and the 
other parameters are: 

(id 



M 



£b [1 + £& 2n~l 

Qo Va 



p(n-2) 
*b 2n-l 



I MX 



(12) 



One can show that R/M changes monotonically along the curve Qo/M = constant. Let 
Qo = M, then (|12"D gives the curve b where Qo = M. Inside (figure |l|, for n > 3) or below 
(figure for n = 3) this curve, Q > M. 

The horizon sizes for the charged black hole are defined as the solutions of equation 

e x (R) = and the results are R± = M ± J M 2 - Q 2 for a given M and Q . Let e A (&) = 0, 



(|) gives curve a where R = R + . Above this curve, we have R < R + , so the solutions 
are unacceptable. We notice this curve is tangential to the curve Qo = M at the point 
£b — ( n — — 1), ct = [(n — l)/(2n — l)] 1_n . On the left hand side of this point, curve a 

represents configurations with their boundary at the external horizon. Consequently physical 
configurations must lie below this curve. 

For given a and one can solve equation @ numerically to get Y(£). We find that the 
central pressure Y(0) will diverge as increases to some value for fixed a < [(n — 1)/ (2n — 
l)] 1_n , then we get the curve c Y(0) = inf which approaches the tangential point of curves 
a and b. Above this curve, there are pressure singularities inside the configurations. So the 
regular static configurations must lie below this curve. (Curves d and e are for later use). 

We see from the figures that, for these special models, no regular solution exists which 
describes a compact charged sphere having Q < M and radius R arbitrarily close to R + 
since any such solution would lie above curve c where a pressure singularity emerges inside 
the sphere. A more general result is assured by the following: 
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Theorem: if the total electric charge of a static perfect fluid ball is smaller than its total 
mass, then there is no regular static configuration having a radius arbitrarily close to the 
size of the external horizon (Yu et al., 1999). 

As stated in (de Felice et al., 1995), the curves in figures |l}|| show that a finite region in the 
parameter space exists, which is the region between curves b and c, where regular solutions 
describing static charged balls of matter may be found with a compactness arbitrarily close 
the critical one for transition to a black-hole state, namely with R —>■ R + and Q — > M. 
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FIG. 4. The pressure of the configurations for n=3 

It is of interest to study the behavior of the pressure inside some of these configurations 
as a function of parameter a. For a fixed n, this parameter gives a measure of the total 
charge within the body, so we show in figure ^ plots of the internal pressure as function of 
radius £, for n = 3, ^ = 0.36 and increasing values of a (the corresponding configurations 
are indicated by '+' in figure g). It is interesting to see how charge contributes to a decrease 
of the central pressure, a property which justifies the very existence of regular solutions 
describing compact balls of matter. 

Regular solutions also exist inside (for n > 3) or below (for n — 3) curve b in figures 1 
and 2, where Qq > M. However, as appears from figure 4, the central pressure is negative 
in these solutions. 

In what follows we shall discuss the stability of these solutions against adiabatic radial 
perturbations. 



III. STABILITY 

We will consider small radial perturbations for the charged, perfect fluid balls, then any 
fluid element which was at r in the unperturbed configuration is displaced to r + £(r, t) in 
the perturbed one where ^<r. In what follows we shall use subscript "i" to denote the 
unperturbed configuration, then we have (Misner et al., 1973): 

r](r,t) = r)i(r) + 5r)(r,t) 
\{r,t) = Ai(r) + 8X(r,t) 
p(r,t) = Pi(r) + 6p(r,t) 
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p(r,t) = pi{r) + Sp(r,t) 
n(r, t) = nj(r) + 5n(r, t) 



(13) 



where n is the baryon number density; the variations are small and such that 5a/ ai ~ £/r, a 
being any of the above parameters. Furthermore we assume that the charge distribution in 
the vibrating configuration remains unchanged, namely Q{r + ^(r, t)) = Qi{r), meaning that 
there are no electric currents for the comoving observer. To first order in the variations, the 
Lagrangian perturbations read: 

Ap(r, t) = p(r + £, t) - pi(r) 

= 6p(r, t) + p^ 
Ap(r,t) = 5p(r,t) + Pi 'Z 
An(r, t) = 5n(r, t) + 
AQ(V,t) = 0. 



From baryon conservation we have: 

dAn 



dr 



—nw 



(14) 



(15) 



Since u r /u l = £, the dot meaning differentiation with respect to the coordinate time t and 
= —1, in the perturbations we have to first order: 

»' = e-Wl-£) 



u 



Using the relation 



■41 



(~g)- 1/2 f(-9) v V 



(16) 



and integrating equation (|T5| ) with respect to coordinate time (first-order analysis), we 
deduce: 



An = —n 
For adiabatic variations we have: 



r- 2 e" A 



' /2 (r 2 en)'4 



(17) 



Ap _ pi 
An n~ 



where 7 is the adiabatic index, then from ( |TTD and (|T^) we get: 



Sp = -j Pi 



r -2 e -A, 



1/2 ( r 2 e A i /2^ + 



From Einstein's equations, we obtain: 



-87r(pi + Pi)e 
7 



-(V + V)e 



(18) 

(19) 
(20) 



thus, integrating fl20D with respect to coordinate time and choosing the integration constant 
so that 5\ = when £ = 0, we get: 



5X = -e Ai 87rr( Pi + Pi )f = -(V + rV)£. 
Using pi]) , (|18D and @ to eliminate respectively 5A, <5p and ^ in (|19D , we have: 



(21) 



(<5p)' = e Ai r 



TT]/ + 1 



<5A 



8^^— - (rV^)' + Stt^V - { Pi + A )]e^ - ^^e 



3 A; 



If we project the identity: 



parallelly to the matter flow, namely: 



T u =0 



u^T\, u = 0, 



(22) 



(23) 



we obtain, from the chosen forms of u M and T^ u : 

dp (p + p) (in Q <iQ 
dr n dr 47rr 4 dr 



0. 



(24) 



Integrating ([24]) with respect to time and using (|14|) and (]17D to remove AQ, Ap and An 
we deduce after some algebra: 



5p = -(pi +pi) 



V A i/2 ( r 2 e V2 e y + sx 



(25) 



If we now project (E3[) transversely to the matter flow, namely: 



where h^ u = 5^ + n^-u^, we obtain: 



c?p 



(p + p)u^u u;fl = -p v - — M y - (T^^.fj) 

CLT 



Oh + Pi )e x ^i = -W - (Sp + s v )\- (a + Pi) 2 47rr4 . 4?rr4 

while the t-component is trivial and gives: 

/ QiQi 



U 



Q dQ 
v 47rr 4 dr 



The r-component of (pq) yields 



(fy)' QiQi' Q| 



/2 



(26) 



e, (27) 



Pi 



47rr4 "(A+ft) T - 



By using the initial value equations fll8|), (|25|) and (p2|) to reexpress 5p, <5p and ((5p)' in 
terms of £ in equation we obtain: 
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(Pi + P l )e^ = { 1Pi r- 2 e*l\' + fo' + ( Pi + + ^"V^C + £( A ' + Pi') 



2 e t/V 



+ (Pi +Pi 



47r 7 p l r-V i+r ' i / 2 C' + ^^e Ai + 47r( Pi + Pi )e Ai £ - 47r Pi 're Ai £ 



QiQi_^, _ Qi f j e »Ji+Ai/2 



4 7rr 4 



4 7rr 4 



7P ^ r- 2 e Ai / 2+3r ' i / 2 C , 



+ e 



Ai/2+3j?i/2 



c 



+ 
+ 



(p? - (Q i QQ/(4 7 rr 4 )) 2 



r2 (Pi + Pi) 



r 3 \ 47rr 4 

QiQi QiQi 



&*Pi(Pi +Pi)^r 7rQi(pi+Pi 



47rr e 



(28) 



where £ = r 2 e is a renormalized displacement function. From the above, the physically 
acceptable solutions of the dynamic equation must satisfy the following boundary conditions: 



i 

r 



finite or zero as 
Ap = - 1Vi r- 2 ^ t2 (rV^y -> 



as 



r -> 
r -> i? 



(29) 
(30) 



where R is the radius of the ball. 

Assume that C( r ^) has a sinusoidal time dependence: 



C(r,t) = C(r)e 



(31) 



hence the dynamic equation and the boundary conditions reduce to an eigenvalue problem 
for the angular frequency uj and amplitude ((r): 



(F(>y + (H + LU 2 W)( = 
so that (/r 3 is finite or zero as r — > and 7j9 i r~ 2 e ,?i//2 C / (r) — > as r — > _R, where: 
F = 1Pl r- 2 e Xi/2+3 ^ 2 



H = e Ai/2+3r?i/2 



A; 



( P ;-(Q i Q0/(47rr 4 )) 2 



r 2 (Pi + Pi) 

Q2f , \ I QiQi QiQi 
i{pi + Pi) + - A — r r 
r u 47rr b 7rr' 



4 / QiQi A s ( 



W = (pi+ Pl )r- 2 e 3 ^ 2+ ^ 2 
By using the variational principle (Mathews J. and Walker R.L., 1970), one gets: 

2 , JoVC 2 " H( 2 )dr 

uj 2 = extreme of J0 v S B 

J R WCdr 

so, the configuration is stable against adiabatic radial perturbations only if: 



(32) 
(33) 

(34) 
(35) 

(36) 
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/ (FC 2 - HC 2 )dr > 
Jo 



(37) 



for any ( which satisfies the boundary conditions fl29|), (p0|). 

Using the simplest trial function ( oc r 3 (Chandrasekhar, 1964) and assuming that the 
adiabatic index 7 was constant throughout the spheres, we calculated the stability of such 
configurations. The domain of stable solutions in the parameter space is shown in figures [l] 
and ^ 

curve d for 7=4, the configurations are unstable above this curve, 
curve e for 7=5/3, the configurations are unstable above this curve, 
and the stability curves approach curve c as the adiabatic index approaches infinity. Evi- 
dently this limit corresponds to perfect rigidity. 

With the inclusion of these curves, figures |1] and ^| allow us to identify a region on the pa- 
rameter space which corresponds to solutions which are physically significant, being regular 
and stable against adiabatic radial perturbations. Since these very compact configurations 
are electrically charged, they are probably not describing astrophysical objects but, perhaps, 
classic models of massive elementary particles. In this case it would be desirable to extend 
our analysis to a more detailed moldelling of the interior, since the adiabatic index depends 
critically on the equation of state (Akmal et al., 1998). 



IV. REGULAR SOLUTIONS AND REGULAR EMBEDDINGS 

Given a Riemannian manifold M with Lorentzian metric g, it may be heuristically helpful 
to find embedding diagrams of any subset of M into Euclidean space. More specifically, one 
has to find a surface in a three-Euclidean space which is isometric to a given two-surface 
of M ( Bini et al., 1999). This problem does not always have a solution. It has been 
argued (Kristiansson et al., 1998) that the limiting condition for regular embedding of space- 
like sections of the vacuum Schwarzschild and Reissner-Nordstrom solutions, conformally 
reduced by the factor Q = g 00 , does coincide with the limiting condition for their regular, 
static, spherically symmetric internal solutions. We shall show here that this is not true. 

Consider the Reissner-Nordstrom space-time solution: 

d # = -(l-^ + ^)d1?+(l-— + Ql) 1 dr 2 + r 2 (d6 2 + sin 2 Qd<\> 2 ) (38) 
then consider the conformal t =constant space-like subspace: 

^=(1-— + ^] ~ dr 2 + (l - — + ®D \ 2 (d6 2 + sin 2 ^ 2 ) (39) 
The two-dimensional surface 9 = n/2, namely: 




2/ "'+gyw do) 



is matched to a surface z = z(R) in the three-dimensional Euclidean space: 
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dR z + R 



2 j i 2 



(41) 



where R, <fi, z are cylindrical coordinates. Comparing (|4l|) with ( |40|) we obtain the condition 
for regular embedding as: 




' i _2m + QT 



1 - 



3m 2Q^ 



- 1 > 



(42) 



From this we deduce, for each M and Qo ; the condition for regular embedding as r > r 
where r is the solution of the equation: 



i _2m + Q| = / i _3m + 2Q|^ 



(43) 




Q /m 

FIG. 5. Configurations in the parameter space for different n and enlargement 

If Qq = 0, we easily deduce f = 9M/4, a well known limit for regular internally static, 
electrically neutral, spherically symmetric and perfect fluid solutions (Buchdahl, 1959). If 
Qo ^ 0, there are several internal solutions which match the external Reissner-Nordstrom 
one. Here, of the various solutions we shall consider those with n = 3 and n = oo. The plots 
of the limit of regularity for these solutions, namely where the central pressure diverges, 
are compared in the (r — Qo)-plane and shown in Figure [| where r is the coordinate radius 
of the sphere's boundary. Although the two curves coincide at r = 9M/4 when Qo = 
and r = M when Qo = M, as expected, they differ from each other everywhere else, curve 
n = oo being at larger radii than curve n = 3 for each value of Qo- Now it happens that 
the condition for infinite central pressure in the n = oo case (a perfect conductor) can be 
expressed analytically through solving the Oppenheimer-Wolkoff equation for conductive 
incompressible perfect fluids and reads: 



\ _ 2m + Ql 



rm 



Q 



3rm - 2QI 



(44) 
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This coincides with condition ( [43|) for regular embedding. It is then clear that regular 
solutions for n = 3 exist which have radii smaller than the corresponding limit of regular 
embedding; this is sufficient to falsify the conjecture. 

V. CONCLUSIONS 

The use of plots in Figures [I] to |5] highlights the result that regular solutions of Einstein's 
equations exist describing charged spheres arbitrarily close to a black-hole. These solutions 
fill a region in a parameter space which is thin but finite so it is not a set of measure zero. 
However not all these solutions are stable against adiabatic radial perturbations. We have 
made a perturbation analysis and have established a stability criterion given by relation (37). 
We have deduced the stability curves for configurations with adiabatic indeces 7 = 5/3 and 
7 = 4 (respectively curves e and d in figures (1) and (2)), assuming that 7 was constant 
throughout the sphere. We see that, for moderate values of 7, the limits of stability prevent 
stable compact charges from being arbitrarily close to a black-hole state but also from having 
a total charge Qo arbitrarily close to their mass. These situations would be possible if the 
adiabatic index diverged; in this case, in fact, the stability curve approaches curve c which, 
we recall, is the limit of regularity for our solutions. Clearly, a diverging 7 corresponds to a 
complete rigidity while, in a less extreme case, 7 would critically depend on the equation of 
state (Akmal et al, 1998). 

Although our solutions can be thought of as being classic (non quantum) models of 
charged particles, it would be interesting to find solutions of Einstein's equations which 
describe charged spheres with a more detailed physical description of their interior. 
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